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Abstract 



1 Introduction 

A selection mechanism in front propagation can be thought of as follows: a certain 
t:J" ' phenomenology is described through an equation that admits an infinite number of 

22 ' traveling- wave solutions, but there is only one which has a physical meaning, the 

one with minimal velocity. Under mild assumptions on initial conditions, the solution 



converges to this minimal- velocity traveling wave. The most remarkable example of this 



ff^ ■ fact is the celebrated F-KPP equation (for Fisher, Kolmogorov-Petrovskii-Piskunov) 

dv 1 d'^v ^, , ^ 

X: ^* ^^^ (1) 

C^ , v{0,x) = Vo{x), X G M. 

Assume for simplicity that / has the form /(s) = s^ — s, but this can be generalized 
up to some extent. We also restrict ourselves to initial data vq that are distribution 
functions of probability measures in M. The equation was introduced in 1937 \^2U\ [2H] 
as a model for the evolution of a genetic trait, and since then has been widely studied 
(in fact more than two thousand works refer to one of these papers). 

Both Fisher and Kolmogorov, Petrovskii and Piskunov proved independently that 
this equation admits an infinite number of traveling wave solutions of the form v{t, x) = 
Wc{x — ct) that travel at velocity c. This fact is somehow unexpected from the modeling 
point of view. In words of Fisher |20l p. 359] 
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Selection Principles. 

"Common sense would, I think, lead us to believe that, though the ve- 
locity of advance might be temporarily enhanced by this method, yet ul- 
timately, the velocity of advance would adjust itself so as to be the same 
irrespective of the initial conditions. If this is so, this equation must omit 
some essential element of the problem, and it is indeed clear that while a 
coefficient of diffusion may represent the biological conditions adequately 
in places where large numbers of individuals of both types are available, it 
cannot do so at the extreme front and back of the advancing wave, where 
the numbers of the mutant and the parent gene respectively are small, and 
where their distribution must be largely sporadic". 

Fisher proposed a way to overcome this difficulty, related to the probabilistic rep- 
resentation given later on by McKean [35], weaving links between solutions to ([1]) and 
Branching Brownian Motion. The general principle behind is that microscopic effects 
should be taken into account to properly describe the physical phenomena. With a 
similar point of view in mind, Brunet, Derrida and coauthors [13^ [14} [TT| [12] started 
in the nineties a study of the effect of microscopic noise in front propagation for equa- 
tion ([T|) and related models, which resulted in a huge number of works that study the 
change in the behavior of the front when microscopic effects are taken into account. 
These works include both numerical and heuristic arguments |13[ [T4l [TTl [T2] [26] as 
well as rigorous proofs [H [5l [161 [311 [32] . Before that, Bramson et.al [TO] gave the ffist 
rigorous proof of a microscopic model for ([T|) that has a unique velocity for every initial 
condition. They also prove that these velocities converge in the macroscopic scale to 
the minimum velocity of ([T|), and call this fact a microscopic selection principle, as 
opposed to the macroscopic selection principle stated above, that holds for solutions 
of the hydrodynamic equation. 

The theory of quasi-stationary distributions (QSD) has their own counterpart. It 
is a typical situation that there is an infinite number of quasi-stationary distributions, 
but the Yaglom limit (the limit of the conditioned evolution of the process started from 
a deterministic initial condition) selects the minimal one, i.e. the one with minimal 
expected time of absorption. 

Up to our knowledge, despite of a shared feeling that similar principles do occur in 
the context of QSD and of traveling waves, this relation has never been stated precisely. 
The purpose of this note is to show that they are two faces of the same coin. We first 
explain this link through the example of Brownian motion. Then we show how to 
extend these results to more general Levy processes. 

The paper is organized as follows. In Section [21 we introduce traveling waves and 
QSDs as macroscopic models. We focus in particular on the KPP equation and the 
links between its traveling waves and the QSDs of a drifting Brownian motion. In 
Section [3l we introduce particle systems enlightening the selection principles observed 
for the macroscopic models. Finally in Section [J] we export these observations to more 
general models. We consider general Levy processes under suitable assumptions and 
analyze selection principles in this context. 

2 Macroscopic models 

We elaborate on the two macroscopic models we study: front propagation and QSD. 
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2.1 Front propagation in the KPP 

Since the seminal papers [201 EH], equation ([T]) has received a huge amount of attention 
for several reasons. Among them, it is one of the simplest models explaining several 
phenomena that are expected to be universal. For instance, it admits a continuum of 
traveling wave solutions that can be parametrized by their velocity c. More precisely, 
for each c G [\/2r, +00) there exists a function Wci M — >• [0, 1] such that 

v{t,x) = Wc{x — ct) 

is a solution to ([1]) . For c < \/2r , there is no traveling wave solution, [H [28] . Hence c* = 
\/2r represents the minimal velocity and Wc* the minimal traveling wave. Moreover, if 
vq verifies for some < & < v2r 






ViTaJ^'il-voix)) = a>0, 
then 



lim v{t, X + ct) = Wc{x), for c = r/b + -b, (2) 

t— ^-oo 2 

see [351 ISS]- If the initial measure has compact support (or fast enough decay at 
infinity), the solution converges to the minimal traveling wave and the domain of 
attraction and velocity of each traveling wave is determined by the tail of the initial 
distribution [H [35l [39] . A smooth traveling wave solution of ([1]) that travels at velocity 
c is a solution to 

-w" + cw + r{w — w) = Q. (3) 

The behavior at infinity of these traveling waves is given by 

^ \c2xe-''^^ c = V2^. 
This behavior is determined by the linearization of ([3]) at tf = 1, i.e. the solution of 

-w" + cw' + rw = 0. (4) 

See |251l39j. We come back to this equation when dealing with quasi-stationary distri- 
butions. 

2.2 Quasi-stationary distributions 

Quasi-stationary distributions have been extensively studied since the pioneering work 
of Kolmogorov (1938), Yaglom (1947) and Sevastyanov (1951) on the behavior of 
Galton- Watson processes. 

The beginning of this theory and an important part of the research in the area 
has been motivated by models on genetics and population biology, where the notion of 
quasi-stationarity is completely natural to describe the behavior of populations that are 
expected to get extinct, conditioned on the event that extinction has not yet occurred, 
on large time scales. 
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Being more precise, consider a Markov process Z = (Zt, t > 0), killed at some state 
or region that we call 0. The absorption time is defined by r = inf{t > 0: Z^ G 0}. 
The conditioned evolution at time t is defined by 

fi]{-):=F^{Zte-\T>t). 

Here 7 denotes the initial distribution of the process. A probability measure u is 
said to be a quasi-stationary distribution (QSD) if /.ij' = v for all t > 0. 

For Markov chains in finite state spaces, the existence and uniqueness of QSDs 
as well as the convergence of the conditioned evolution to this unique QSD for every 
initial measure follows from Perron-Frobenius theory. The situation is more delicate 
for unbounded spaces as there can be 0, 1 or an infinite number of QSD. Among those 
distributions, the minimal QSD is the one that minimizes E^(t). 

The Yaglom limit is a probability measure v defined by 

u := lim /ij", 

t— >oo 

if it exists and does not depend on x. It is known that if the Yaglom limit exists, then 
it is a QSD. A general principle is that the Yaglom limit selects the minimal QSD, 
i.e. the Yaglom limit is the QSD with minimal mean absorption time. This fact has 
been proved for a wide class of processes that include birth and death process, Galton- 
Watson processes, random walks, Brownian motion, more general Levy processes, etc. 
It can also be proved for i?-positive processes by means of the theory of i2-positive 
matrices [38]. To give a flavor of the results that hold in this situation, consider a 
discrete time Markov chain in N that it is absorbed at 0. Denote p = {p{i,j), i,j G N) 
its transition matrix so that p is sub-stochastic. We use p^"^ for the n-th power of 
p. We say that p is i?-positive if one (and hence both) of the following equivalent 
statements hold 

1. For some i and j, the sequence R^p^'^>{i,j) tends to a finite non-zero limit as 

n — )• 00. 

2. There exist non-negative, non-zero eigenvectors u = {l'{k))k^^., /3 = (/3(A;))fcgN 
associated to the eigenvalue 1/R such that X^fc^i ^{^)l^{^) < 00. 

In 1966, Seneta and Vere-Jones proved the following theorem 

Theorem 2.1 (Seneta and Vere-Jones, |38j). Assum.e that the m.atrix p is R-positive, 
then the conditioned evolution converges to u as n ^ 00 if one of the following condi- 
tions hold 

1. The left eigenvector v satisfies ^i^{i) < 00 and the initial distribution /_f is dom- 
inated (pointwise) by a multiple ofu. 

2. The right eigenvector (3 is bounded away from zero and ^ /u(j)/3(j) < 00 

Observe that in both situations we have that u is the minimal QSD and the Yaglom 
limit. Also every initial distribution with tail light enough is in the domain of attraction 
of v. So, in the i?-positive case, the situation is pretty clear. These results can be 
applied for instance to the Galton- Watson process. In that case, a detailed study of 
the domain of attraction of the other QSDs (which are parametrized by an interval) is 
given in [37] where it can be seen that the limiting conditional distribution is given by 
the tail of the initial distribution. 
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Unfortunately, on the one hand i?-positivity is a property difficult to check and 
on the other hand, there is a lot of interesting processes that are not i?-positive. 
For example, a birth and death process with constant drift towards the origin, has 
a continuum of QSDs and initial distributions with light tails are attracted by the 
minimal QSD, which can be computed explicitly [19], but this process is not i?-positive 
and hence Theorem 12.11 does not apply. 

The presence of an infinite number of quasi-stationary distributions is something 
anomalous from the modeling point of view, in the sense that no physical nor biological 
meaning has been attributed to them. The reason for their presence here and in the 
front propagation context is similar: when studying for instance population or genes 
dynamics through the conditioned evolution of a Markov process, we are implicitly 
considering an infinite population and microscopic effects are lost. 

So, as Fisher suggests, in order to avoid the undesirable infinite number of QSD, 
we should take into account microscopic effects. A natural way to do this is by means 
of interacting particle systems. We discuss this in Section [3l 

Brownian Motion with drift. Quasi-stationary ditributions for Brownian Mo- 
tion with constant drift towards the origin are studied in [33^ ^^ . We briefly review 
here some of the results of these papers and refer to them for the details. 

For c > we consider a one-dimensional Brownian Motion X = {Xt)t>o with drift 
— c defined by Xt = Bt — ct. Here Bt is a one dimensional Wiener process defined in the 
standard Wiener space. We use ¥x for the probability defined in this space such that 
Bt is Brownian Motion started at x and Ej, for expectation respect to P^- Define the 
hitting time of zero, when the process is started at x > by Tx{c) = inf{t > 0: Xt = 0} 
and denote with Pt the submarkovian semigroup defined by 

Pi7(^)=IEx(/(^t)l{.4c)>t})- (5) 

In this case, differentiating ([5]) and after some manipulation it can be seen that the 
conditioned evolution fiPt has a density u{t, •) for every t > and verifies 

du 1 3 VL dii 1 3vb 

u{t,0) = n(i,+oo) = 0, t > 0, 

Recall now that a probability measure u in M+ is a QSD if 

F,{Xt G-\Xt>0) = !/(•). 

It is easy to check that if i/ is a QSD, the hitting time of zero, started with v is an 
exponential variable of parameter r and hence z^ is a QSD if and only if there exists 
r > such that 

uPt = e~''V, for any t > 0. 

Differentiating ([5]) and using the semigroup property we get that i/ is a QSD if and 
only if 

j\\f"-cf')du = -r j fdu, forall/GCo^(M+). (7) 
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Integrating by parts we get that the density w of i' must verify 

-w" + cw' + rw = 0. (8) 

Solutions to this equation with initial condition w{0) = are given by 

{me~^^ sin(-v/c^ — 2rx) r > ^, 
me~'^^ sinh(vc?--2rx) r <^. 

Observe that w defines an integrable density function if and only if < r < c^/2 (or 
equivalently, c > v2r ) • One can thus parametrize the set of QSDs by their eigenvalues 
r, {vr : < r < c^/2}. For each r, the distribution function of u^., v{x) = f^ w{y) dy is 
a monotone solution of ([S]) with boundary conditions 

u(0) = 0, ^(+oo) = 1, (9) 

which is the same equation Q but in a different domain. The following theorem 
characterizes the domain of attraction of each QSD. 

Theorem 2.2 (Martinez, Picco, San Martin, [33]). Let j be a probability measure on 
(0, +oo) with density p. If 

lim \ogp{x) = b < c, (10) 

x/"oo X 

then liuit pj = i^^ffe); where r{b) = cb — 5^/2. 

Observe that this last equation is equivalent to c = r/b + ^b, which should be 
compared with ([2]). Also remark that if (jlOp holds then 



lim log/_i([x, +oo) = 6. 

x/'oo X 

We come back to equation ([8]) later, shedding light on the links between QSD and 
traveling waves. 

Finally let us mention that a similar result holds for the discrete-space analog, i.e. 
birth and death processes with constant drift towards the origin [171 HI] 



3 Particle systems 



In this section we introduce two particle systems. The first one is known as Branching 
Brownian Motion (BBM) with selection of the A*" right-most particles (A— BBM). As a 
consequence of the link between BBM and F-KPP that we describe below, this process 
can be thought of as a microscopic version of F-KPP. The second one is called Fleming- 
Viot and was introduced by Burdzy, Ingemar, Holyst and March [15], in the context 
of Brownian Motion in a d-dimensional bounded domain. It is a slight variation of 
the original one introduced by Fleming and Viot [21] . The first interpretation of this 
process as a microscopic version of a conditioned evolution is due to Ferrari and Marie 

m. 
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3.1 BBM and F-KPP equation 

One-dimensional supercritical Branching Brownian Motion is a well- understood object. 
Particles diffuse following standard Brownian Motion started at the origin and branch 
at rate 1 according to an offspring distribution that we assume for simplicity to be 82- 
When a particle branches, it has two children and then dies. As already underlined, its 
connection with the F-KPP equation and traveling waves was pointed out by McKean 
in the seminal paper [35] • Denote with Nt the number of particles alive at time t > 
and ^i(l) < • • • < Ct(-^f ) the position of the particles enumerated from left to right. 
McKean's representation formula states that if < ^0(2;) < 1 and we start the process 
with one particle at (i.e. A^(0) = 1, ^o(l) = 0), then 

vit,x):=E(^voiUi)+x)] 

is the solution of ([T]). Of special interest is the case where the initial condition is the 
Heaviside function vq = 1{[0, +00)} since in this case 

vit,x) = p(ei(i) + X > 0) = p(et(iVt) < x). 

This identity as well as various martingales obtained as functionals of this process have 
been widely exploited to obtain the precise behavior of solutions of ([1]), using analytic 
as well as probabilistic tools P [9] [2il [25l [35| [39] . 

3.2 A/^— BBM and Durrett-Remenik equation 

Consider now a variant of BBM where the A^ right-most particles are selected. In 
other words, each time a particle branches, the left-most one is killed, keeping the 
total number of particles constant. 

This process was introduced by Brunet and Derrida |1H [T2] as part of a family of 
models of branching-selection particle systems to study the effect of microscopic noise 
in front propagation. By means of numerical simulations and heuristic arguments, they 
conjectured that the linear speed of A— BBM differs from the speed of standard BBM 
by (log A^)~^ and in a series of papers with coauthors they study various statistics of 
the process |13[ [T^ [TT] I12j . Berard and Gouere proved this shift in the velocity for a 
similar process. We refer to the work of Maillard [31] for the most detailed study of 
this process. 

Durret and Remenik [16] considered a slightly different process in the class of Brunet 
and Derrida: A— BRW. The system starts with A particles. Each particle gives rise 
to a child at rate one. The position of the child of a particle at x G M is x + y, where 
y is chosen according to a probability distribution with density p, which is assumed 
symmetric and with finite expectation. After each birth, the A + 1 particles are sorted 
and the left-most one is deleted, in order to keep always A particles. They prove that 
the empirical measure of this system converges to a deterministic probability measure 
Vf for every t, which is absolutely continuous with density u(t, •), a solution of the 
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following free-boundary problem 
Find {■~f,u)such that 

du 



oo 



^^(*'^) = / uit,y)pix-y)dy Vx > 7(t), 



u{t,y)dy = 1, u{t,x)=0, Vx < 7(t), 

7(t) 

u{0,x) = uo{x). 

They also find all the traveling wave solutions for this equation. Just as for the BBM, 
there exists a minimal velocity c* G M such that for c > c* there is a unique traveling 
wave solution with speed c and no traveling wave solution with speed c for c < c* . The 
value c* and the behavior at infinity of the traveling waves can be computed explicitly 
in terms of the Laplace transform of the random walk. In Section U] we show that these 
traveling waves correspond to QSDs of drifted random walks. 

It follows from renewal arguments that for each A^, the process seen from the left- 
most particle is ergodic, which in turn implies the existence of a velocity vn at which 
the empirical measure travels for each N. Durrett and Remenik prove that these 
velocities are increasing and converge to c* as A^ goes to infinity. 

We can interpret this fact as a weak selection principle: the microscopic system has 
a unique velocity for each A^ (as opposed to the limiting equation) and the velocities 
converge to the minimal velocity of the macroscopic equation. The word "weak" here 
refers to the fact that only convergence of the velocities is proved, but not convergence 
of the empirical measures in equilibrium. 

In view of these results, the same theorem is expected to hold for a A^— BBM that 
branches at rate r. In this case the limiting equation is conjectured to be given by 

Find {'j,u)such that 

du 1 d u 

-qI^^^^) = 252;^*^*'^^ + ™^*'^^ Vx>7(t), ^^2) 

u{t,y)dy = 1, n(t,x)=0, Vx < 7(t), 



00 



7(i) 

m(0, x) = uq{x). 

The empirical measures in equilibrium are also expected to converge to the minimal 
traveling wave. More precisely, 

Conjecture 3.1. Both N-BBM and N-BRW are ergodic, with (unique) invariant mea- 
sure X^ and the empirical measure distributed according to X^ converges to the delta 
measure supported on the minimal quasi- stationary distribution. 

Traveling waves. Let us look at the traveling wave solutions u{t,x) = w{x — ct) 
of (112p . Plugging-in in (I12p we see that they must verify 

(13) 

Which is exactly ([8]) . Note nevertheless that in (|13p the parameter r is part of the data 
of the problem (the branching rate) and c is part of the unknown (the velocity) , while 



1 // / f°° 

-w" + cw -|- rtf = 0, w(^) =0, / w(y) dy = 1. 

2 Jo 
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in ([8]) the situation is reversed: c is data (the drift) and r unknown (the absorption 
rate under the QSD). However, we have the following relation 



c is a minimal velocity for r r is a maximal absorption rate for c 

in ([13D ^^ in (IHD 



Observe also that 1/r is the mean absorption time for the QSD associated to r 
and hence, if r is maximal, the associated QSD is minimal. So the minimal QSD for 
Brownian Motion in M+ and the minimal velocity traveling wave of (|12p are one and 
the same. They are given by 

Uc*(r)ix) = Ur*(c) = 2r*xe~^^'"* = (c*)^xe~'^ ^, 

which is the one with fastest decay at infinity. 

Again, the distribution function t" of u is a monotone solution to the same problem 
but with boundary conditions given by v{0) = 0, t;(+oo) = 1. 

It is worth noting that although the solutions to ^ and (j3|) are not the same since 
they are defined in different domains, there is a natural way to identify them (and also 
with solutions of ([3])). Given positive constants c and r, there is a solution u; of (l3|) if 
and only if there is solution w oi ^. Moreover, we have 

hm ^ = 1. 

x-5>oo w{x) 

In this case, the proof of this statement is immediate since solutions of ([8]) and ^ are 
explicit and the relation among solutions of (j3]) and ([3]) is very well understood |25| . 
We conjecture that the same situation holds in much more generality. 

3.3 Fleming- Viot and QSD 

The Fleming- Viot process can be thought of as a microscopic version of conditioned 
evolutions. Its dynamics are built with a continuous time Markov process Z = (Zf, t > 
0) taking values in the metric space Au{0}, that we call the driving process. We assume 
that is absorbing in the sense that 

¥{Zt = 0\Zo = 0) = 1, Vt > 0. 

We use r for the absorption time 

T = inf{t >0: Zt^ A}. 

As, before, we use Pt for the submarkovian semigroup defined by 

Ptf{x) = EM{Zt)l{T>t}). 

For a given N > 2, the Fleming- Viot process is an interacting particle system with N 
particles. We use ^t = (6(1)) • • • )Ct(^)) ^ ^^ to denote the state of the process, ^(i) 
denotes the position of particle i at time t. Each particle evolves according to Z and 
independently of the others unless it hits 0, at which time, it chooses one of the A^ — 1 
particles in A uniformly and takes its position. The guenine definition of this process is 
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not obvious and in fact is not true in general. It can be easily constructed for processes 
with bounded jumps to 0, but is much more delicate for diffusions in bounded domains 
[6l [22] and it does not hold for diffusions with a strong drift close to the boundary of 

A, 13. 

Here we are also interested in the empirical measure of the process 



N 



1 ^ 

^E%W- (14) 



1=1 

Its evolution is mimicking the conditioned evolution: the mass lost from A, is redis- 
tributed in A proportionally to the mass at each state. Hence, as A'" goes to infinity, we 
expect to have a deterministic limit given by the conditioned evolution of the driving 
process Z, i.e. 

^if (A) ^ F{Zt € A\t > t) (N^oo). 



This is proved in [30] by the Martingale method in great generality. See also [23] for 
a proof based on sub and super-solutions and correlations inequalities. A much more 
subtle question is the ergodicity of the process for fixed N and the behavior of these 
invariant measures as A^ — )■ oo. As a general principle it is expected that 

Conjecture 3.2. If the driving process Z has a Yaglom limit u, then the Fleming-Viot 
process driven by Z is ergodic, with (unique) invariant measure X^ and the empirical 
measures (|14p distributed according to X^ converge to v. 



We refer to [23] for an extended discussion on this issue. This conjecture has been 
proved for subcritical Galton- Watson processes, where a continuum of QSDs arises [2]. 

We have again here a microscopic selection principle: whereas there exists an in- 
finite number of QSDs, when microscopic effects are taken into account (through the 
dynamics of the Fleming-Viot process) , there is a unique stationary distribution for the 
empirical measure, which selects asymptotically the minimal QSD of the macroscopic 
model. 

When the driving process is a one dimensional Brownian motion with drift — c 
towards the origin as in Section 12. 2( the proof of the whole picture remains open, but 
the ergodicity of FV for fixed A^ has been recently proved [3l [27] . 

So, from [401 Theorem 2.1] we have that for every t > 0, /if^ converges as A^ — )• oo 
to a measure fXt with density u{t,-) satisfying ([6|). The open problem is to prove a 
similar statement in equilibrium. Observe that u is a stationary solution of ([6]) if and 
only if it solves ([8|) for some r > 0. Hence, although equations (fT2]) and Q are pretty 
different, stationary solutions to ([6|) coincide with traveling waves of (112p . 

3.4 Summing up 

1. The link between A^— BBM and Fleming Viot, in the Brownian Motion case is 
clear. Both processes evolve according to A^ independent Brownian Motions and 
branch into two particles. At branching times, the left-most particle is eliminated 
(selection) to keep the population size constant. The difference is that while 
A^— BBM branches at a constant rate Nr, Fleming-Viot branches each time a 
particle hits 0. This explains why in the limiting equation for A^— BBM the 
branching rate is data and the velocity is determined by the system while in the 
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hydrodynamic equation for Fleming- Viot the velocity is data and the branching 
rate is determined by the system. 

2. The empirical measure of A^— BBM is expected to converge in finite time intervals 
to the solution of (jl2|) . This is supported by the results of [TB] where the same 
result is proved for random walks. 

3. The empirical measure of Fleming- Viot driven by Brownian Motion converges in 
finite time intervals to the solution of ([6]). 

4. Both A^-BBM seen from the left-most particle and FV are ergodic and their 
empirical measure in equilibrium is expected to converge to the deterministic 
measure given by the minimal solution of (I13p . 

Note though that while for iV— BBM r is data and minimality refers to c, for 
Fleming- Viot c is data and minimality refers to 1/r (microscopic selection prin- 
ciple). 

5. u{t, x) = w{x — ct) is a traveling wave solution of ()12p if and only if w is the 
density of a QSD for Brownian Motion with drift — c and eigenvalue — r. 

6. c is minimal for r (in (I13p ) if and only if 1/r is minimal for c. So, we can talk of a 
"minimal solution of ()13p ". which is both a minimal QSD and a minimal velocity 
traveling wave. 

7. The microscopic selection principle is conjectured to hold in both cases, with the 
same limit, but a proof is still unavailable. 

4 Traveling waves and QSD for Levy processes 

Let Z = {Zt, t > 0) be a Levy process with values in M, defined on a filtered space 
(ri, J-", (J-i),P) and Laplace exponent -0 : M — )■ M defined by 

E(e^^') = e'^(^)* 

such that 

where 6 G M, o" > (which ensures that Z is non-lattice) and g is defined in terms of 
the jump measure 11 supported in M \ {0} by 

g{9) = f{e'^ - 1 - fel{|,|<i|)n(dx), / (1 A x2)n(dx) < oo. 

Jx JR 

Let 9* = sup{0: |'(/'(^)| < oo} and recall that ip is strictly convex on (0,0*) and by 
monotonicity 'il){9*) = 'ip{6*—) and tp'{6*) = %Ij'{6*—) are well defined as well as the 
right derivative at zero ^''(O) = ^"'(0+) = lE(Zi), that we assume to be zero. We 
also assume that 9* > 0. In this case we can relate ^ to the characteristic function 
^(A) = - logE(e'^^i) by i^{9) = -^{-19) for < 6* < 6*. 

This centered Levy process plays the role of Brownian Motion in the previous 
sections. 

The generator of Z applied to a function f ^ Cq, the class of compactly supported 
functions with continuous second derivatives, gives 

Cf{x) = la^f"{x) + bfix) + / (fix + y)- f{x) - ?//'(x)l{|y| < l})U{dy). 
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The adjoint of C is also well defined in Cg and has the form 

C*f{x) = la^f"{x) - bf'ix) + / (fix -y)- f{x) + ?//'(x)l{|y| < l})U{dy). 

Now, for c > we consider the drifted process Z^ given by 

Zf = Zt-ct 

It is immediate to see that the Laplace exponent of Z^ is given by V'c(^) = ^(^) — c9 
for 9 £ [0, 9*], that Cq is contained in the domain of the generator Cc of Z'^, and that 
Ccf = Cf — cf. Recall that the forward Kolmogorov equation for Z is given by 

j^E^f{Zt)) = Cf{x), 

while the forward Kolmogorov (or Fokker-Plank) equation for the density u (which 
exists since cr > 0) is given by 

— u(t,x) = C*u{t,-){x). 

As in the Brownian case, we consider 

• A branching Levy process (BLP) (Nt, (^t(l)i • • • ,S,t{J^t))) driven by C. 

• A branching Levy process with selection of the N rightmost particles (A^— BLP), 
also driven by C. 

• A Fleming- Viot process driven by Cc (FV). 

We focus on the last two processes. For a detailed account on BLP, we refer to [29]. 
Let us just mention that the KPP equation can be generalized in this context to: 

dv 

— = Cv + rf{v), xGM, t>0, 

^* (15) 

v{0, x) = Vo{x), X G R. 

A characterization of the traveling waves as well as sufficient conditions of existence 
are then provided in |29j . 

For A^— BLP we expect (but a proof is lacking) that the empirical measure converges 
to a deterministic measure whose density is the solution of the generalized Durrett- 
Remenik equation 



C*u{t,x) + ru{t,x), X > 7(t), ,,„x 

1, u{t,x) = 0, X < 7(t), 

uo{x), x >0. 

Existence and uniqueness of solutions to this problem have to be examined, but we 
believe that the proof presented in [16] can be extended to this context. 



Find 


{-l,u) 


such that 




/■OO 


du 




/ 


u{t,y)dy 
u{0,x) 
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We show below the existence of travehng wave solutions for this equation under 
mild conditions on C based on the existence of QSDs. 

Concerning FV, it is known [50] that the empirical measure converges to the deter- 
ministic process given by the conditioned evolution of the process, which has a density 
for all times and verifies 

— (i,x) = C*u{t,x) + c—{t,x)-u{t,x)C*u{t,y,)dy t>0,x>0, 
u(t,0) = u{t,+oo) = 0, t > 0, 

4.1 Traveling waves vs. QSDs 

We now establish the link between traveling waves and QSD. 

Proposition 4.1. The following statements are equivalent: 

• The probability measure v with density w is a QSD for Z'^ with eigenvalue —r , 

• n(t, x) = w{x — ct) is a traveling wave solution with speed c for the free-boundary 
problem (J16p . with parameter r. 

Proof. Denote {f,g) = J f{x)g{x)dx. A QSD v for Z'^ with eigenvalue — r is a solution 
of the equation 

{yC, + rv,f)=Q,^feCl 

Using that £*/ = C* f + cf and writing that i/ has density w, we obtain that 

£*w + cw' + rvu = 0, 

which in turn is clearly equivalent to w being a traveling wave solution with speed c 
for {[161). ■ 



4.2 Minimal QSD and minimal velocities 

We now study the relation between minimal traveling waves and minimal QSD. Using 
known results on QSDs for Levy processes, we can describe the traveling waves for the 
corresponding set of equations. Since we rely on results of Kyprianou and Palmowski 
|30j . we assume that Z is non-lattice and consider the following two classes. 



Definition 4.1. The Levy process Z belongs to class CI if there exists < Oq < 9* 
such that iP'{Oq) > and the process (Z, P^°) is in the domain of attraction of a stable 
law with index 1 < a <2, where the probability P^ is given by 



dP^ 

dp: 



^e{Zt-z)-i>(e)t _ 



Tt 



Definition 4.2. The process Z is in class C2 if — oo < if)' {9*) < 0, and the function 
X —7- IIq*[[x,oo)) is regularly varying at infinity with index —f3 < —2, where Ilg{dx) = 
e^^n(dx). 
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Recall that we are also assuming V^'(O) = and hence there exists a critical c*, 
possibly infinity, such that for c < c* , ipc attains its negative infimum at a point 9c < 9* 
with ipc{9c) = and for c > c* , the negative infimum is attained at 9*. In this case we 
write 9c = 9\ Since ^^{0) = '4'{9) - c9, observe that e^'^^t-^)-Md)t ^ ^e{Zt-z)-^{9)t _ 
Hence, if Z is in class CI, then for c < c* , Z*^ is in class A in the sense of Bertoin and 
Doney |12| [30] . Similarly if Z is in class C2 and c > c* , then Z^ is in class B in the 
sense of Bertoin and Doney. 

The union of class A and B represents a very large family of Levy processes including 
for instance Brownian motion and spectrally negative (positive) processes, as well as 
many others. The following theorem is proved in |30| . 



Theorem 4.1 (Kyprianou and Palmowski, [30]). Assume c < c* and Z is in class CI 
or c> c* and Z is in class C2. Then the Yaglom limit of Z'^ exists and is given by 

u{dx) = 9QKeS^,9c)e-'^^''VeSx)dx. 

Here kq and Vq are respectively the Laplace exponent of the ascending ladder process 
and the renewal function of the ladder heights process corresponding to {Z'^,¥^). 

The proof is based on a careful control of the asymptotics of the process as i — )• oo 
that in particular yields 

P^(r > t) ~ H{x, 9c)£{t)e^^''^^^^ . (18) 

Here H is a function that depends on the characteristic exponent of the process and i 
is a function regularly varying at infinity. This allows us to prove the following 



Proposition 4.2. The probability measure v defined in Theorem, \J^.1\ is the minimal 
QSD. 

Proof. Recall that under a QSD, the hitting time of zero is exponentially distributed. 
Since ly is the Yaglom limit, for any bounded function / 



/ 



/d.= li,n'^'<^™-^>'> 



t^oo P^(r > t) 
Applying this to f{y) = Py(T > s) and using Corollary 4 in [30], one obtains that 

P.(r > .) = lim ^i^f'l^"^ = eMM0c)s). 

t->oo Fa;(T > t) 

Hence the parameter of v equals —tpd^^c)- H ^ is not minimal, there exists another 
QSD u, with parameter f > —ipc{9c)- Then Ej>(e^'^'''^''-''^) < oo and as a consequence 
there exists x > such that 'E^{e~'^''^^''^^) < oo, but this contradicts ([T8|). 



Let C be the subset of M^ such that (c, r) G C if and only if there exists a QSD v 
for Z'^ with eigenvalue — r. Proposition 14.11 states that this set coincides with the set 
of pairs (c, r) such that there exists a traveling wave for (|16p with velocity c. 

Definition 4.3. We say that r is maximal for c if r = max{r': {c,r') £ C}. In the 
same way, c is minimal for r if c = min{c': {c',r) € C}. 



Selection Principles. 



15 



Proposition 4.3. Under the same hypotheses of Theorem ]^. l\ we have that c is minimal 
for r if and only if r is maximal for c. 

Proof. For each c > 0, the niaxhxial absorption rate is given by r{c) = —ipdOc)- The 
proof follows by observing that the function c i— )■ '(/'c(^c) is strictly decreasing and 
continuous. I 



5 Conclusions and general conjectures 

We emphasized the direct links between quasi stationary distributions QSDs for space 
invariant Markov processes in M and traveling waves of Durrett-Remenik equation, 
that are closely related to the generalized F-KPP equation. We proved also that 
minimal QSDs correspond to minimal velocity traveling waves. As a general fact both 
macroscopic and microscopic selection principles are expected to hold for QSDs and 
for traveling waves. They have been proved for a series of models. This suggests that 
the selection principles in front propagation and in QSDs are one and the same. For 
random walks in M, the microscopic selection principle is an open problem in both 
cases (traveling waves and QSD). 
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